THE RIEMANN PROBLEM FOR THE STOCHASTICALLY 
PERTURBED NON- VISCOUS BURGERS EQUATION AND THE 
PRESSURELESS GAS DYNAMICS MODEL 



ANASTASIYA KORSHUNOVA, OLGA ROZANOVA 

Abstract. Proceeding from the method of stochastic perturbation of a Langevin 
system associated with the non-viscous Burgers equation we construct a so- 
lution to the Riemann problem for the non-interacting particles and sticky 
particles systems. We analyze the difference in the behavior of discontinuous 
solutions for these two models and relations between them. 



Let us consider the Cauchy problem for the non- viscous Burgers equation: 

dtu + {u,\7)u — 0, u{x,0) = uq{x), (1) 

where u{x,t) — {ui, is a vector-function R"+-'^ W\ 

It is well known that on the smooth solutions of the Burgers equation this equa- 
tion is equivalent to the system of ODEs 

x{t) ^u{t,x{t)), u{t,x{t)) = (2) 

for the characteristics x = x(t). 

We associate with ([2]) the following system of stochastic differential equations: 

dXk{t) = Uk{t)dt + nd{Wk)u 
dUk{t)^Q, fc = l,...,n, (3) 
X{0) = X, U{0) = u, 
where {X{t),U{t)) runs in the phase space M" x M", tr > is a constant and 
{W)t — {Wk)ti k = 1, n, is the n - dimensional Brownian motion. 
Let us introduce a function 

/ uP{t, X, u)du 
J F{t,x,u)du 

where P{t, x, u) is the probability density in position and velocity space. This value 
has a sense of the conditional expectation of U for fixed position X Ij . If we choose 

n 

Pq{x,u) = S{u - uo{x))fo{x) = Jl S{uk - {uo{x))k)fo{x), (5) 
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where foix) is an arbitrary sufficiently regular nonnegative function sucli that 
/ fo{x)dx = 1, then u{0,x) = uo{x). Certain properties of u{t,x) was studied 

R" 

in [2] and [3]. 

The density P = P{t, x, u) obeys the Fokker-Planck equation 



dP 
'dt 



V- d A 1 2 52 

. k=l k=l 



(6) 



subject to initial data ([5]). 

We apply the Fourier transform in the variables x and w to ([6]) and ([5]) and 
obtain the Cauchy problem for P = Pit, A, ^): 



which solution is given by the following formula: 

P(t,A,O-mA,^ + A0e-^^'l^l''. (9) 
The inverse Fourier transform allows to find the density P{t,x,u): 

P{t,x,u) = J J e^(^'-)e'(«'«)pdAde = 

R" R" 

h.0(3)t+3-x|^ 

S[u — uo[s)) jo{s) e 2^^t ds. (10) 



K" 

Then we substitute P{t, x, u) in ^ and get the following expression for ■u{t, x): 

J ■"o(s)/o(s)e ds 

u{t, x) = — ^ . (11) 

J Jo[s)e ds 

R" 

It should be noted that integrals in pTjl are defined also for a wider class of /o (x) 
then the probability density of the particle positions in the space at the initial 
moment of time. If the integral / fo{x)dx diverges (for example, for fo{x) = const), 

R" 

we consider the domain [—L, L]", where L > and use another definition of u{t, x): 

/ Mo(s)/o(s)e ds 

u{t,x)= Mm tMl^ . (12) 

J /o(s)e as 

Evidently, this definition coincides with (fTTj) for fo{x) e ii(M"). 
The following property of ■u{t, x) holds: 

Proposition 1. Let uo{x) and fo{x) > be bounded functions of class and 
the solution to the respective Cauchy problem {!]) keeps smoothness for t < t^ < 
+00. Then u(t,x) tends to solution of problem as a —>■ for any fixed {t,x) e 
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Proof. Let us denote J{uq{x)) the Jacobian matrix of the map x i — > uo{x). As it 
was shown in |4j (Theorem 1), if J{uq{x)) has at least one eigenvalue negative for 
a certain point x £ R", then the classical solution to ([T]) fails to exist beyond a 
positive time t*. Otherwise, t^. = oo. The matrix C{t,x) = (/ + tJ{u{){x))), where 
/ is the identity matrix, fails to be invertible for t — t^. 
The formula ((Tl|) (or ([12])) imphes 

/ Wo(s)/o(s) Inn ^-^^L-^e'^^^^'B^ ds 

\im u{t,x) = ^ |„„(,),+,_,|2 • 

r /o(s) lim ^ , . e ds 

If p{t, X, s) = ito(s)t + s — X , we can use locally the implicit function theorem 
and find s = s(t, x,p). Moreover, dp — det C{t, s) ds. Therefore, 

/ Ms{p))fo{s{p)) detiC{t,s{p)))-'Sip)dp 

<T^o } fo{s{p)) Aet[C{t,s{p))) ^5{p)dp 

where we denote by so(i, x) — s{t, x, 0) the vector-function which obeys the follow- 
ing vectorial equation: 

UQ{so{t,x))t + so{t,x) — X = Q. (13) 
Let us show that u{t,x) — UQ{so{t, x)) satisfies the Burgers equation, that is 

n n 

^ 5j(uo,i)(soj)t + ^ uoj'9fe(uo,i)so,/c^^ = 0, z = l,...,n, (14) 

and uo(so(0, x)) — uq{x). Here we denote by uo,i and so,i the i - th components of 
vectors mq and sq, respectively. 

We differentiate (fT5)) with respect to t and xj to get the matrix equations: 

n 

^ (soj)f -I- uo,j = 0, i = l,...,n, 

and 

n 

Qfc (scfc)^;^ + (^ij = 0, i,j = 1, ...,n, 

fc=i 

where (5^ is the Kronecker symbol. The equations imply 

n 

{so,j)t = - ^ (C"^).y uo,j, {sQ,k)x, = -{C^^)jk, j = l,..,n. (15) 

i=l 

It remains now only to substitute ^T5\i into p^ . 
Further, implies uo(so(0, 2:)) = uq{x). □ 

It is important to note that so(t,x) is unique for all t for which the solution to 
the Burgers equation u{t, x) is smooth. 
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Let US denote p{t,x) — J P{t,x,u)du. From (fTO|) we have 

R" 

p{t,x) = , ^ , / /o(s)e S^'^ ds. (16) 



(V27rtcr)' 

It can be readily checked that in the one-dimensional case the functions p{t, x) 
and x) solve the following system; 

dtp + d^ipu) = icT^ d.^^p, (17) 
5t(pM) + 9,(pu2) ^ _ _ ufd^Pdu + iff^ a^^(pt,). (18) 



The equation (|T7|) follows from the Fokker-Planck equation (|6]) directly. To prove 
(|17p we note that definitions of u(t, x) and /3(i, x) imply 

dt{pu) — dt / uP{t,x,u)du= / uPl.{t,x,u)du 



u'^P^{t,x,u)du + —cr'^ I uP'J.^{t,x,u)du 

u^P^it, X, u)du + icr2a^^(pu). 
Further, 

dx{pu^) — u^dxP + 2pudxU 
= 2udx{pu) — u^dxP ~ J 2uuP!^{t,x,u)du^ J u^P^{t, x,u)du. 

R R 

A certain analog of system (fT71) . (ll8p was obtained in [S]. 

Let us denote f{t,x) = lim p{t,x). Taking into account Proposition 1 and the 

Fokker-Plank equations, as a limit cr — > for smooth f{t,x) and u(t,x) we obtain 
the system of pressureless gas dynamic (e.g. [5]) in any space dimensions: 

dtf + divxifu) = 0, 

(19) 

dt{fu) + V{fu<»u)=0. 

However the formula pip has sense also for discontinuous initial data (/o (a;), uo{x)). 
For the sake of simplicity we dwell on the one-dimensional case and consider the 
following initial data: 

foix)=fi+f20{x-xo), (20) 

Uq{x) = Ui + U20{x — Xq), (21) 

where 0{x ~ xq) is the Heavisidc function. Without loss of generality we assume 
Xq = 0. 
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Definition 1. We call the couple of functions {f(t,x),u{t,x)) the generalized so- 
lution to the problem 119\) I120\) . \21\) in the sense of free particles, if for almost all 
(i, x) e R+ X M" 

f{t,x) — lim(lini u{t,x) — lim{lmi (t, x)), 

where (p^{t,x),u^{t,x)) satisfy the system |J7| j, US]) with initial data (/q (a;), tig (a;)) 
from the class C^(]R") such that lim/Q(x) = fa{x), limuQ(x) = uq(x) for almost 

all fixed X e M." . 

Proposition 2. The solution {f{t,x),u{t,x)) to the problem U9\) . i20\} . i21\) in the 
sense of free particles does not depend of the choice of the couple (/q(x), 7/5(0;)) G 
Ci(R")nC6(R"). 

Proof. Let us choose two couples of smoothed functions {fQi{x) , Uqi{x)) and 
(/o2(a;),"o2(2;)) such that 

hni /o^i (x) = hni f^2 (x) = fo (x) , hm ul^ix) ^ Una (x) = uo{x) 

£—►0 e-^O e^O e^O 

for any fixed x S M". Then the couple 

if^ix),uUx)) = (/o^(x) - /o^2(x), u'olix) - "02(2;)) e Ci(R) n a(K") 

can be considered as initial data for the problem (ll9p - (f2T|) . To prove the proposition 
we have to show that the respective solution is identically zero. 
Indeed, from (fTH)) we have 



f{t,x) = lim / f^{s)S{s - sn{t,x))ds = \im{f^,{so{t,x)) - f^^{so{t, x))) = 0. 

Here so{t,x) is a solution of (fTH]) as before. Analogously proceeding from pT|) . we 
prove that u{t,x) = 0. □ 

Our purpose is to find relations between a stable solution to the Riemann problem 
of system that can be obtained as a limit ct ^ from (|17|1. (fT8|) (as we will show 
below it is not necessarily looks like (fT9l) !) and the couple {f{t,x),u{t,x)). 

The Riemann problem in ID case. According to Definition 1 we must con- 
sider the smoothed initial data instead of (pO|) and pTjl . As follows from Proposition 
1 we can choose any couple of smoothed initial data we want. It will be convenient 
to consider the piecewise linear approximation of initial data of the form 



foix) = ii^x + fi + i^, -e<x<e, (22) 




"0(2^) = \ TT^ + "1 + TT' —£< X < e, (23) 



where /i, /2, wi and U2 are constants. 

Note that these functions can be pointwisely approximated by functions from 
the class C^(W). 
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From (fTH]) we can find the density p'^{t, x) corresponding to the smoothed initial 
data (/o^(x),M§(x)): 



1 



(24) 



where $(a) = — / e "2 dx is the Gauss function, Ci = uit — x — e, 
V27r _oo 



(iti + U2)i — a; + £, and 

ir = F^{t,x) 



o 



(25) 



The expression for I^''^ can be written out, however, we are interested only in 
behavior of I^'"^ as cr ^ 0. It can be calculated that lim F^{t,x) — 0. 



To find u{t,x) we compute the numerator in formula 

Wo(s)/o(s)e ^^^^ ds 



/27rtcr 



Uip''(t,x) +M2(/l +/2)«' 



^2 ) 



where 



ir^N^t,x) 



O cr^ + ere" 



(26) 



where N^it-, x) 



{x ~ (lii + — )t) + — 1 F'(t, x) and lim 7V^(t, x) = 0. 

U2t + 2e ^ 2 ' ' 2 I ^ ' e^a ^ ' 



Thus, we have the following result: 



it^(i, x) = ui + 



e.(T 
2 



(27) 



where ' and 1 2 are given by in (|25p and (pS)) . respectively. Note that > ±00 

fT 

as CT ^ 0. 

Now we can find the generalized solution to the Riemann problem as 
f{t, x) — lim ( lim Pe{t^ x)), u(t, x) — lim ( lim Ue{t^ x)). 

Let us introduce the points x\ = uit — e and xl = (ui +U2)t + e. Their velocities 
are ui and ui + M2, respectively. 
We consider two cases: 

1. M2 > (velocity of the point x^ is higher than velocity of the point x\). We 
can find f^{t,x) — lim p^{t,x) from Let us note that this formula contains 

F'^{t,x). It is easy to see that 

lim x^ — X — uit, lim xf, — x — (ui + U2)t 
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and \im F^t,x) = 0. Thus, 





f{t,x) = lim r{t,x) 



/l, 

2 ' 
0, 

/1 + /2 
2 ' 
/1 + /2, 



re < Ult, 

X = Ult, 

Ult < X < {ui + U2)t, 
X = {Ui+ U2)t, 
X> {Ui+ U2)t. 



initial data u^{t, x) = lim u^{t, x) as follows: 



Further, from ([27l) we find the solution of the gas dynamic system with smooth 



X X ^ '-^2' 
X ^ X^ . 



Ul, 



u'^{t, x) 



F{ey 
U1+U2, 



Ul 



It can be shown that 

Nie) 



lim 



F{e) 



lim ( — 

e— 



"2 



U2t + 2e 



[x 



(-1 + f )t) 



X 

--Ul. 



Thus, we get the following solution: 



Ml, 

X 



u(t,x) — Mmu'^ {t, x) = \—, 



X < Ult, 

Ult < X < (ui + U2)t, 



(28) 



Ul + U2, X > {Ui + U2)t. 



We can see that the velocity includes the rarefaction wave. It is well known stable 
solution to the Burgers equation with initial data pT|) . 

It is interesting to note that if we compute the limit in e first we get the solu- 

U2 

tion u(t,x) = Ul + U20{x — {ui H — which is unstable with respect to small 
perturbations. 

2. U2 < (the velocity of il is higher than the velocity of if). From 
(P7)l we find as before: 



and 



f{t,x) 



fi, 

3/1 + /2 

2 ' 

2/1 + /2, 
3/1 + 2/2 

2 

/1+/2, 



a; < (ui + U2)t, 

X ^ {Ui+ U2)t, 

{ui + U2)t < X < Ult, 

X ~ Ult, 

X > Ult, 



Ul, 



u{t, x) = lui + 



fi + h 



X < {ui+ U2)t, 
-U2, {ui + U2)t < X < Ult, 
X > Ult. 



2/1 + /2 
^Ui + U2, 

Remark We can consider in this framework the singular Riemann problem with 
initial density fo{x) = /i + f20{x - xq) + fsSix). 
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The Hugoniot conditions and the spurious pressure. As we have been 
proved, if / and u are smooth, they solve the pressureless gas dynamics system. 
Now we ask the question which system satisfy the solution of this system with 
jumps in the sense of free particles. 

The system of conservation laws ([19]) implies two Hugoniot conditions that should 
be held on the jumps of the solution This signifies the solution satisfies the 
system in the sense of integral identities. If we denote by D the velocity of the 
jump and — h{y + 0) — h[y — 0) the value of the jump, then the continuity 

equation and the momentum conservation give [f]D = [fu] and [fu]D = [fu^], 
respectively. 

In the case U2 > the velocity is continuous, therefore the Hugoniot conditions 
hold trivially. 

We should check these conditions for the jumps in the case U2 < 0. An easy 
computation shows that the first one is satisfied, however, the second one does not 
hold. To understand the reason let us estimate the integral term in (fT8|) on the 
generalized solutions in the case U2 < : 

{u — u{t, x))'^ Px{t, X, u) du — 



^ ^ fo[s){uo{s)-u{t,x)f ( e-'""*t."^r"' ) da 



R 

2 / _ (uo(s)t + s-a;)2 

/o(s)((uo(s)-WFp(t, a:))+(uFp(i, a;)-u(i,a;))) e ds = 



R 

-2{uFp{t,x) - u{t,x)) J_ [ fois){uois)-UFp{t,x)) (e"' L-^t -\ ds+ 



'2'Kta 

R 

{uFp{t,x) ~ u[t,x)) — / /o(sj e ds 



'l-Kta 

R 

The integrals I2 and /a tend to zero as cr ^ due to properties of the Riemann 
data since u(t^ x) —> UFp(t, x) for almost all a: € R. Let us estimate Ii. 



h ^--jL^ J fi{ui-u{t,x)f (^e-'"'2t^r' ^ ds- 

{h + f2){{ui+U2)-u{t,x) f ( e-""'*"'^^^""' ) ds = 







1 +/2)^M2 ( -("1*-^)" 



/2^a iVi + h? 
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Thus, 

{^ji + h) 

Thus, the integral term corresponds to a spurious pressure between the jumps 
X = + U2) t and X — uit, namely, 

p = ^j^'/Jff ieix - (^.1 + u,) t) - 9{x - u, t)). (29) 

The Hugoniot condition [fu]D ~ [fu^ + p] is satisfied with this kind of pressure. 
Thus, we get the following theorem. 

Theorem 1. The generalized solution to the Riemann problem i20\) . for 
the pressureless gas dynamics system in the sense of free particles (according to 
Definition 1) in the case of a discontinuous velocity (u2 solves in fact the gas 

dynamics system with a pressure defined by I129\) . 



Sticky particles model vs non-interacting particles model. In our model 
the particles are allowed to go through the discontinuity as one particle does not 
feel the others. However, if we are in the frame of the sticky particles model we 
should assume that the particles meeting one other stick together on the jump. 
The non-interacting particles model and the sticky particles model are equivalent 
for smooth velocities, however, if the velocity has a jump, the solutions behavior 
differs drastically. Nevertheless, we can study the solution to the Riemann problem 
in the case of U2 < for the sticky particles model, too, basing on the solution 
obtained in the present paper. Indeed, the jump position x{t) is a point between 
= (f^i + ""2) t and X2it) = ui t. The mass m{t) accumulates in the jump due 
to the impenetrability of the discontinuity with the velocity 

m{t) = {x{t) - {ui + U2)t)ii2f, + /2) - /i) + (u.t ~ x{t))ii2fi + f^) - (/i + /a)) 

= + U2){fl + /2) - fl)t + X{t)f2) = -[Uf]t + [f]x{t), 

where [ ] stands for a jump value. Further, if we change heuristically the overlapped 
mass between xi and X2 to the mass concentrated at a point, then from the condi- 
tion of equality of momenta in the both cases we can find the velocity of the point 
singularity: 

(Ml + W2)((2/l + /2) - fl){x{t) - {Ul + U2)t) + Ul fl {Ult - Xit)) 

^-[u^f]t+[uf]xit)^mit)xit). 
Thus, to find the position of the point singularity we get the equation 

{[f]x{t) - [uf]t) ±{t) = [uf]x{t) - [u^f]t, (30) 
subject to initial data a;(0) = 0. The respective solution is 

x{t) = [[uf] ± v/["/]^-[/]N2/]) t, (31) 

where the sign should be chosen from the condition xi(t) < x{t) < X2(t). It can be 
readily shown that the latter condition is satisfied either for plus or minus in the 
formula (|3ip . The condition coincides with the Lax stability condition ui < x{t) < 

Ul + U2- 
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The formulas describing the amplitude of the delta-function in the density com- 
ponent and the singularity position obtained earlier in [8].[9].[l0] give the same 
result. 

It is worth mentioning that the spuriovis pressure (|29p does not arise in the sticky 
particles model. 

Back to the Burgers equation. Now we are able to get the solution to the 
Cauchy problem (|2ip . Let us assume that fo{x) = const, that is /2 = 0. Then 
for the case U2 < from ([50]) we obtain the well known formula for a position of 
the jump: 
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If U2 > 0, the solution is continuous and it is given by the formula (pS)) . 
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